Idempotents, Mattson-Solomon Polynomials and Binary 

LDPC codes 

R. Horan, C. Tjhai, M. Tomlinson, M. Ambroze and M. Ahmed 
Fixed and Mobile Communications Research 
University of Plymouth, PL4 8AA, United Kingdom 
February 1, 2008 

Abstract 

We show how to construct an algorithm to search for binary idempotents which 
may be used to construct binary LDPC codes. The algorithm, which allows control 
of the key properties of sparseness, code rate and minimum distance, is constructed 
in the Mattson-Solomon domain. Some of the new codes, found by using this 
technique, are displayed. 



1 Introduction and Background 

The use of idempotents in the construction of cyclic error correcting codes is well es- 
tablished and the resulting literature is extensive (for example, see [1], [3], [7]). The 
basic building blocks for this theory are the primitive idempotents. Any cyclic code 
may be described by a unique idempotent and this idempotent is a sum of primitive 
idempotents. For binary cyclic codes, efficient algorithms exist for the calculation of 
these primitive idempotents. 

Another way of constructing idempotents in the binary case is by using cyclotomic 
cosets and it was this property which was exploited by Shibuya and Sakaniwa in [4]. 
Their goal was to use idempotents to construct parity check matrices for LDPC codes 
which have no cycles of length 4 in their factor graphs. At the heart of their technique 
is a lemma which is a variation of a result used by Weldon [8], for the construction of 
difference set cyclic codes. Using this lemma and a subsequent theorem, they were able 
to simplify the problem of determining which of the idempotents that are constructed, 
using a single cyclotomic coset, do not have cycles of length 4. They then extended this 
theory to more general idempotents. 

This approach to the construction of LDPC codes has the great advantage of simplic- 
ity, the parity check matrices depend only upon the correct choice of cyclotomic cosets 
and these are very easily calculated. However, we believe that this advantage is offset 
by some fundamental weaknesses. 

Whilst the absence of 4 cycles is a desirable objective in the construction of LDPC 
codes it is not mandatory [5], since there are some good codes which do not have this 
property. An example of such a code is included in this paper. The code rate is also an 
important property of codes but, as Shibuya and Sakaniwa admit in their conclusion, 
the codes which they construct in this way are 11 expected to have a large minimum 
distance at the expense of rate 11 (our italics). The minimum distance of a code is a 
crucial property but there is no indication in [4] of how either a single cyclotomic coset, 
or combinations of more than one cyclotomic cosets, should be chosen to guarantee that 
the code constructed has a large minimum distance. 



In order to address the question of how to choose idempotents which will produce 
good LDPC codes we propose an entirely different route. As in [4], we shall deal ex- 
clusively with binary cyclic codes. Making effective use of the Mattson-Solomon poly- 
nomial, we produce an algorithm which not only allows us to choose, in a systematic 
way, idempotents with low weight, and therefore a correspondingly sparse parity check 
matrix, but also with the desirable features that the corresponding codes have a high 
code rate and a large minimum distance. 

This paper is organised as follows. In section 2 we shall review the necessary the- 
ory and explain how it will be used to provide an algorithm for the determination of 
idempotents which may be used to construct good codes. In section 3 the design and 
implementation of this algorithm is given and then, in section 4, some of the results are 
displayed. Finally, in section 5, we draw our conclusions on this approach. 

2 Binary Idempotents 

Let F = GF(2), n be a positive integer and T be the splitting field for x n — 1 over F. 
Let a G T be a primitive nth root of unity and let T(x) be the polynomials in J-[x] of 
degree < n — 1. If a(x) G T(x) then the map $ : T — > T is defined by 

Ma)](z)=J2^> n - J (1) 

and $(a) is the Mattson-Solomon polynomial of a (see [1]). (We use x and z for the poly- 
nomial variables to distinguish between the polynomials in the domain and codomain 
of $.) If o is multiplication of polynomials mod (x n — 1) and * is defined on T(z) by 
the rule (^aiZ 1 ) * (J^biZ 1 ) = ^a^z 1 then it is well known [1],[7], that 

$ :(T,+,o)^(T,+,*) 

is an isomorphism of rings, in particular it is an isomorphism of the additive groups. 

If S(x) is the subset of T(x) consisting of polynomials with coefficients in GF(2) 
(binary polynomials) and E(x) is the subset of T(x) consisting of idempotents, both of 
these subsets are additive subgroups of T(x). It is easy to show (see [1]) that 

*:(£(*),+) - (E(z),+) (2) 
$ :(E(x),+) - (S(z),+) (3) 

are both isomorphisms and from this it is obvious that 

$ : (S(x) n E(x), +) - (E(z) n S(z), +) (4) 

is also an isomorphism. 

Suppose that u(x) is a binary idempotent which is used to construct a parity check 
matrix for a cyclic code. The parity check matrix is constructed from the n-cyclic shifts 
of u(x) [2], and so for the resulting code to be a LDPC code, u(x) must have low weight. 

If h(x) = gcd(a: n — l,u(x)) and g(x) = (x n — \)/h(x), then g(x) is the generator 
of the cyclic code. If the generator, g(x), has degree n — k, the dimension of the code 
is k and the larger the value of k, the better the code rate. Since g(x) is a divisor of 
x n — 1, all of the zeros of g(x) are nth roots of unity, and there are n — k of these. 



Further, gcd(g(x), h(x))—l and x n — 1 = h(x)g(x), so that the number of distinct nth 
roots of unity which are also roots of u(x) is k. The dimension of the code is therefore 
the number of nth roots of unity which are also roots of u(x). 

The BCH bound of the code is determined by the number of consecutive powers of 
a, taken cyclically (modn), which are also roots of g(x). For the reasons outlined in 
the previous paragraph, this is precisely the same as the number of consecutive powers 
of a, taken cyclically (modn), which are not roots of u(x). 

The important features of the code are therefore determined by: 

(a) the weight of the idempotent u(x), 

(b) the number of nth roots of unity which are roots of u(x), 

(c) the number of consecutive powers of a which are not roots of u(x). 

Take u(x) € S(x) fl E(x) and let = 9 be its MS polynomial. The inverse 

mapping 

S" 1 : (S(z) n E(z), +) -> (E(z) n S(x), +) (5) 

is defined as follows: If A(z) = [<J>(a)](z) is the Mattson- Solomon polynomial of the 
polynomial a(x) = a + a±x + • • • + a n x n ~ x then, for i — 0, . . . , n — 1, 

a % = -A(a l ) (6) 
n 

(see [1]). Let h(z) = gcd{9(z),z n - 1) and let f(z) = {z n - l)/h(z). The three key 
properties relating to the idempotent u(x), listed above, are easily gleaned from its 
Mattson-Solomon polynomial 9{z), and f(z), as follows: 



2.1 The weight of u(x) 

The weight of u(x) is the number of nth roots of unity which are zeros of f(z). To see 
this note that f(a l ) = if and only if 9(a l ) = 1, since idempotents take only the values 
and 1 in T . Now u = and the coefficients of u(x) = Uo + U\x + . . . + Un-ix 11 ' 1 

are given by 

u i = 6{a i ) mod 2 for i = 0, . . . n - 1 (7) 

(cf equation (6)). Thus Ui = 1 precisely when f(a l ) = 0, giving the weight of u(x) as 
the degree of the polynomial f(z). 

2.2 The zeros of u(x). 

From the definition of the MS polynomial, (1), 

n 

e( z ) = u(a j )z n ~ j (8) 
j'=i 



and the number of zeros of u(x) which are roots of unity is clearly n — wt(9(z)). 



2.3 The BCH bound of the code. 



The BCH bound of the code is the largest number of consecutive powers of a which are 
not roots of u(x), i.e. the number of consecutive i, taken (modn), such that u{a l ) = 1. 
From (8), this is the largest number of consecutive non-zero coefficients in 9, taken 
cyclically (mod n). 

Using this information, a systematic search for idempotents can now be made in increas- 
ing order of weight, with accompanying knowledge of the number of roots which are nth 
roots of unity and the corresponding BCH bound. This algorithm is constructed in the 
Mattson-Solomon domain. 

Let the decomposition of z n — 1 into irreducible (over F = GF(2)) polynomials be 
z n - 1 = h{z)f 2 {z) . . . f t {z). For i = 1, . . . , t, let h{z) = (z n - l)/f t (z) and let 9 t (z) 
be the associated primitive idempotent (see [1] or [7]). These are displayed below in an 
array, together with other idempotents: 

ui(x) 
u 2 (x) 

u t (x) 

Here u\(x) , U2{x) , . . . , u t (x) are the idempotents whose Mattson-Solomon polynomials 
are 9i(z), 02(z), . . . , 9 t (z), respectively. Let / C {1, 2, . . . , t} and let u, 9 and / be defined 
as u = J2 ieI Ui, 9 = J2 ieI 9i and f{z) = Yl ieI fi(z). From the properties of primitive 
idempotents, if h(z) = gcd(9(z), z n — 1) then it follows that gcd(f(z),h(z)) = 1 and 
z n — 1 = f(z)h(z). The idempotent u will now have the following properties. 

wt(«) = ^deg(/ 4 ), (10) 

number of zeros of u = n — wt{9) . (11) 

The BCH bound is determined from 9(z) as explained in 2.3. 

Since methods for finding the 9i and fa are well documented (see e.g. [6]) a search 
algorithm can be built around this observation to find a suitable weight idempotent 
with a known number of zeros and a known BCH bound. The rows of the array (9), 
are ordered by the degree of the polynomials, i.e. deg(/j) < deg(/ i+ i) for all i, and a 
search can be made in increasing order of weight. When a successful outcome has been 
obtained, only at this stage is the inverse Fourier transform (MS -1 ) evaluated to find 
the corresponding idempotent. All of the information which is required will already be 
known. 

3 Design and Implementation 

If t denotes the number of cyclotomic cosets modulo n, the complexity of an exhaustive 
search algorithm is 0(2*). We reduce this search complexity by targeting the search on 
the three key parameters: 



6i(z) fi(z) 

&2{Z) f 2 (z) 

0t(z) f t (z) t 



3.1 Sparseness of the parity-check matrix 



In [8], Weldon introduced difference-set cyclic codes. These codes have the desirable 
property that the parity check equations are orthogonal on all bits and have no cycles 
of length 4 in their factor graphs. A necessary condition for this is that if v(x) is the 
polynomial which generates the parity check matrix then the weight of v(x) must satisfy 
the inequality 

wt(v(x))(wt(v(x)) - 1) < n, (12) 

where n is the code length. Since the weights of the idempotents u(x) are related to the 
degrees of the /j by (10), a reasonable bound is 

£deg(/i)<v^. (13) 

iCI 

In practice we have gone a little beyond this limit and this has enabled us to find some 
good codes which do have cycles of length 4 in their factor graph. 

3.2 Code-rate 

The code-rate is directly proportional to the number of roots of u(x). If we let R m in 
represent our minimum desired code-rate then, following equation (11), we can refine 
our search bound to 

wt(0) < (1-R min )n. (14) 

3.3 Minimum distance 

Let d be the lowest desired minimum distance and let rg be the largest number of 
consecutive non-zero coefficients, taken cyclically mod n, of 6. Then, following the 
discussion of 2.3, we restrict our search algorithm to those 9 for which 

r e > d (15) 

We develop an efficient, but exhaustive recursive tree-search based on the above bounds. 
The developed search algorithm, Algorithm 1, is initialised by setting V and index to 
and —1 respectively. 

4 Code Example and Performance 

Since the algorithm is an exhaustive search, the code construction method presented in 
this paper is able to produce, in addition to new codes, many well-know cyclic codes, 
for example the Difference-Set Cyclic codes and the Euclidean and Projective Geometry 
codes. Some of the new codes which we have found using this technique are presented in 
Table 1. All codes in Table 1, except those labelled with *, have orthogonal parity-check 
equations. 



Algorithm 1 CodeSearch(V, index, F(x), Q(z)) 



Input: R m i n <= minimum code-rate of interest 
d <= lowest expected minimum distance 
5 <= small positive integer 

F(x) -<= {fi(x)} Vi G / sorted in ascending order of the degree 
Q(z) <= {Oiiz)} Mi el 
Output: CodesList contains set of codes 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 



T <= V 

for (i—index+1; i < Size(J); i++) do 
T -i= T 

1 prev 1 

if {J2\/jeT de M( x )) + de S(fi( x )) <Vn + $) then 
Append i to T 

Ev,- 6 t^(*) 
if (wt(9(z)) < (1 — R min )n and r e > d) then 
u(x) <= MS' 1 (9{z)) 
if u(x) is non-degenerate then 
C -<= a cyclic code defined by 
if (C £ CodeList) then 

Add C to CodeList 
end if 
end if 
end if 

CodeSearch(T, index, F(x), Q(z)) 
end if 

J- ^ J- prcv 

end for 



(n, k) 


u{x) 


drain 


(51,26)* 


1 +x 3 +x 6 +x 12 +x 17 +x 24 +x 27 +x 34 +x 39 +x 45 +x 48 


10 


(63,44)* 


1 + X 7 + X 9 + X 14 + x 18 + x 27 + x 28 + X 35 + X 36 + 
2^45 _|_ ^49 _|_ ^,54 _|_ ^,56 


8 


(93,47) 


1 + x 2 + X 8 + X A1 + X 62 + x 35 + x 47 


8 


(105,53) 


1 + x 4 + x 3U + x 32 + x 45 + x 46 + X 53 


8 


(117,72)* 


1 + X + x 2 + x 4 + x 8 + x 11 + x 16 + x 22 + x 32 + x 44 + 

^,59 _|_ ^64 _|_ ^»88 


12 


(127,84)* 


1 + X + x 2 + x 4 + X 8 + x 16 + x 32 + X 55 + X 59 + x 64 + 
x 91 _|_ ^,93 _|_ ^,109 _|_ ^,110 _|_ ^,118 


10 


(219,101) 


1 + x 2 + x 8 + x 32 + x Y3 + x 74 + x n + x 89 + x iiU + 
x 128 + X 137 


12 


(255,135) 


1 + x 4 + x 13 + x 21 + X 39 + x 54 + X 55 + X 91 + X 121 + 
X 123 _|_ x 148 + x 195 


13 


(255,175) 


1 + X + X 3 + X 7 + x 15 + x 26 + x 31 + X 53 + X 63 + X 98 + 
X 107 -4- x 127 + x 140 + x 176 + x 197 + x 215 


17 


(341,205) 


1 + x 29 + X 87 + x 92 + x 94 + x 114 + x 122 + x 156 + x 202 + 

x 203 _|_ ^,213 _|_ ^,217 _|_ ^,234 _|_ ^,257 _|_ ^,273 


16 


(511,199) 


l + x + x 3 +x 7 + x 15 +x 31 +x 63 + x 82 + x 1U0 + x 127 + 
x 152 -4- x 165 + x 201 + x 255 + x 296 + x 305 + x 331 + x 403 


19 


(511,259) 


1 + x 31 + x ±2 + x 93 + x 115 + x 217 + ^240 + ^261 + 

-360 _i_ ™420 _i_ ™450 _i_ ™465 
X ■+- X ■+- X ■+- X 


13 


(819,435) 


1 + X + X 3 + X 7 + x 15 + x 31 + x 63 + x 127 + x 204 + 

x 255 _|_ ^409 _|_ x 5ll 


13 


(819,447) 


1 + X + X 3 + X 7 + X 15 + X 31 + X 63 + X 127 + X 204 + 

x 255 + x 350 + x 409 + x 511 + x 584 + x 701 


16 



Table 1: Examples of the constructed codes 
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Figure 1: Frame error performance of the (127,84) cyclic code 



Throughout the paper, it is assume that the codewords are transmitted across a noisy 
communication channel with BPSK modulation and at the receiving end is the modified 
Belief-Propagation decoder which approximates the Maximum-Likelihood decoder [9]. 



Figure 1 shows the frame-error-rate (FER) performance of the (127, 84) cyclic code, 
which is a code which does have cycles of length 4. Neverthless, the performance of this 
code is outstanding and, at 10 -3 FER, it is within 0.2dB of the sphere-packing-bound 
constraint for binary transmission. 
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Figure 2: Frame error performance of the (255, 175) codes 

Figure 2 shows the performance of two (255, 175) codes, one constructed using our 
method and one an irregular computer generated code. We can see that our code, 
which achieves a coding gain of around 0.4dB compared to the equivalent irregular 
code, performs approximately withinr 0.15dB of the sphere-packing-bound constraint 
for binary transmission. 
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Figure 3: Frame error performance of the (341,205) codes 

Our construction method can produce LDPC codes with high minimum distance and 
therefore they do not suffer from error-floor. Figure 3 demonstrates the performance of 
our (341,205) code, which is inferior to the equivalent irregular code in the low signal- 
to-ratio region, but the irregular code exhibits early error-floor due to its low minimum 
distance. 



5 Conclusions 



A method of constructing binary cyclic codes from the finite-field transform (Mattson- 
Solomon) domain is able to produce a large number of codes which have high minimum- 
distance and code-rate. These codes have sparse parity-check matrix and thus are 
applicable as LDPC codes. Due to their cyclic property these LDPC codes have n 
parity-check equations instead oin — k equations as in the case of random LDPC codes. 
With these extra parity-check equations to iterate with, the performance of the iterative 
decoder is improved. 

In designing cyclic LDPC codes of length n, our new method allows one to increase 
the d m i n of the code by combining additional irreducible factors of 1 + z n which, in turn, 
reduces the sparseness of the parity-check matrix. The ability to control the sparseness 
of the parity-check matrix is a trade-off against the minimum distance of the code. 

Simulation results have shown that the our cyclic codes have outstanding perfor- 
mance which are superior to the equivalent irregular LDPC codes. The high d m i n of 
these codes ensures the absence of an early error-floor in their performance. 
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